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Abstract. The deletion-contraction algorithm is perhaps the most popular method for 
computing a host of fundamental graph invariants such as the chromatic, flow, and reli- 
ability polynomials in graph theory, the Jones polynomial of an alternating link in knot 
theory, and the partition functions of the models of Ising, Potts, and Fortuin-Kasteleyn 
in statistical physics. Prior to this work, deletion-contraction was also the fastest known 
general-purpose algorithm for these invariants, running in time roughly proportional to 
the number of spanning trees in the input graph. 

Here, we give a substantially faster algorithm that computes the Tutte polynomial — and 
hence, all the aforementioned invariants and more — of an arbitrary graph in time within a 
polynomial factor of the number of connected vertex sets. The algorithm actually evaluates 
a multivariate generalization of the Tutte polynomial by making use of an identity due to 
Fortuin and Kasteleyn. We also provide a polynomial-space variant of the algorithm and 
give an analogous result for Chung and Graham's cover polynomial. 

An implementation of the algorithm outperforms deletion-contraction also in practice. 

1. Introduction 

Tutte's motivation for studying what he called the "dichromatic polynomial" was algo- 
rithmic. By his own entertaining account [JT], he was intrigued by the variety of graph 
invariants that could be computed with the deletion-contraction algorithm, and "playing" 
with it he discovered a bivariate polynomial that we can define as 

(1) T G (x,y)= J2(x-l) c{F) - ciE \y-l) ciF)+lFHVl . 

FCE 

Here, G is a graph with vertex set V and edge set E; by c{F) we denote the number of 
connected components in the graph with vertex set V and edge set F. Later, Oxley and 
Welsh [36] showed in their celebrated Recipe Theorem that, in a very strong sense, the 
Tutte polynomial Tq is indeed the most general graph invariant that can be computed using 
deletion-contraction. 

Since the 1980s it has become clear that this construction has deep connections to many 
fields outside of computer science and algebraic graph theory. It appears in various guises 
and specialisations in enumerative combinatorics, statistical physics, knot theory and net- 
work theory. It subsumes the chromatic, flow, and reliability polynomials, the Jones polyno- 
mial of an alternating link, and, perhaps most importantly, the models of Ising, Potts, and 
Fortuin-Kasteleyn, which appear in tens of thousands of research papers. A number of sur- 
veys written for various audiences present and explain these specialisations |39[ H3] W$\ 05] . 

*Lund University, Department of Computer Science, P.O.Box 118, SE-22100 Lund, Sweden. E-mail: 
andreas .bjorklundOlogipard. com, thore .husf eldtOcs .lu.se. 

^Helsinki Institute for Information Technology HUT, Department of Computer Science, University of 
Helsinki, P.O.Box 68, FI-00014 University of Helsinki, Finland. E-mail: petteri.kaski@cs.helsinki.fi, 
mikko.koivisto@cs.helsinki.fi. This research was supported in part by the Academy of Finland, Grants 
117499 (P.K.) and 109101 (M.K.). 

1 



2 



Computing the Tutte polynomial has been a very fruitful topic in theoretical computer 
science, resulting in seminal work on the computational complexity of counting, several 
algorithmic breakthroughs both classical and quantum, and whole research programmes 
devoted to the existence and nonexistence of approximation algorithms. Its specialisation 
to graph colouring has been one of the main benchmarks of progress in exact algorithms. 

The deletion-contraction algorithm computes Tq for a connected G in time within a 
polynomial factor of t(G), the number of spanning trees of the graph, and no essentially 
faster algorithm was known. In this paper we show that the Tutte polynomial — and hence, 
by virtue of the Recipe Theorem, every graph invariant admitting a deletion-contraction 
recursion — can be computed in time within a polynomial factor of o~(G), the number of 
vertex subsets that induce a connected subgraph. Especially, the algorithm runs in time 
exp (O(n)), that is, in "vertex-exponential" time, while r(G) typically is exp(u;(n)) and 
can be as large as n n ~ 2 [12J. Previously, vertex-exponential running time bounds were 
known only for evaluations of Tq in special regions of the Tutte plane (x,y), such as for 
the chromatic polynomial and (using exponential space) the reliability polynomial, or only 
for special classes of graphs such as planar graphs or bounded-degree graphs. We provide 
a more detailed overview of such prior work in ^2} 

1.1. Result and consequences. By "computing the Tutte polynomial" we mean com- 
puting the coefficients of the monomials x l y 3 in Tc{x,y) for a graph G given as input. 
Of course, the coefficients also enable the efficient evaluation of Tq(x, y) at any given point 
(x,y). Our main result is as follows. 

Theorem 1. The Tutte polynomial of an n-vertex graph G can be computed 

(a) in time and space o~(G)n°^' ; 

(b) in time 3 n n ^ and polynomial space; and 

(c) in time 3 n ~ s 2 s n°^ and space 2 s n°^ for any integer s, < s < n. 

Especially, the Tutte polynomial can be evaluated everywhere in vertex-exponential time. 
In some sense, this is both surprising and optimal, a claim that we solidify under the 
Exponential Time Hypothesis in §2.5| Moreover, even for those curves and points of the 
Tutte plane where a vertex-exponential time algorithm was known before, our algorithm 
improves or at least matches their performance, with only a few exceptions (see Figure [TJ. 

For bounded-degree graphs G, the deletion-contraction algorithm itself runs in vertex- 
exponential time because t(G) = exp (O(n)). Theorem 111 still gives a better bound because 
it is known that o~(G) = 0{{2 — e) n ) for bounded degree jTj Lemma 6], while t{G) grows 
faster than 2.3™ already for 3-regular graphs (see |2.4[ ). The precise bound is as follows: 

Corollary 2. The Tutte polynomial of an n-vertex graph with maximum vertex degree A 
can be computed in time ^n°^ l \ where £a = (2 A+1 — 1) 1 /( A + 1 ). 

The question about solving deletion-contraction based algorithmic problems in vertex- 
exponential time makes sense in directed graphs as well. Here, the most successful attempt 
to define an analogue of the Tutte polynomial is Chung and Graham's cover polynomial, 
which satisfies directed analogues to the deletion-contraction operations [13 . It turns out 
that a directed variant of our main theorem can be established using recent techniques that 
are by now well understood, we include the precise statement and proof in Appendix [Cj 

1.2. Overview of techniques. The Tutte polynomial is, in essence, a sum over connected 
spanning subgraphs. Managing this connectedness property introduces a computational 
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Figure 1. An atlas of the Tutte plane (x,y). The five points shown by 
circles and the points on the hyperbola (x — l)(y — 1) = 1 are in P, all 
other points are ^P-complete. Those points and lines where algorithms with 
complexity exp(0(n)) were previously known (sometimes only in exponential 
space), are labelled with their running time; note that the hyperbolas (x — 
l)(y— 1) = q were known to be vertex-exponential only for fixed integer q. See 
§2.3| for references. Our result is that the entire plane admits algorithms with 
running time 2 n n°^ and exponential space, or time 3 n n°^ and polynomial 
space. The only points that are known to admit algorithms with better 
bounds are the "colouring" points (—2,0) and (—3,0), the "Ising" hyperbola 



(x — l)(y — 1) = 2, for which a faster algorithm in observed in {2.3 and of 



course the points in P. (Only the positive branches of the hyperbolas are 
drawn.) 



challenge not present with its specialisations, e.g., with the chromatic polynomial. Neither 
the dynamic programming algorithm across vertex subsets by Lawler [31] nor the recent 
inclusion-exclusion algorithm [8], which apply for counting fc-colourings, seems to work 
directly for the Tutte polynomial. Perhaps surprisingly, they do work for the cover polyno- 
mial, even though the application is quite involved; the details are in Appendix [C] and can 
be seen as an attempt to explain just how far these concepts get us. 

For the Tutte polynomial, we take a detour via the Potts model. The idea is to eval- 
uate the partition function of the q-states Potts model at suitable points using inclusion- 
exclusion, which then, by a neat identity due to Fortuin and Kasteleyn |16[ [39] . enables 
the evaluation of the Tutte polynomial at any given point by polynomial interpolation. Fi- 
nally, another round of polynomial interpolation yields the desired coefficients of the Tutte 
polynomial. Each step can be implemented using only polynomial space. Moreover, the 
approach readily extends to the multivariate Tutte polynomial of Sokal [39 which allows 
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the incorporation of arbitrary edge weights; that generalisation can be communicated quite 
concisely using the involved high-level framework, which we do in £[3j To finally arrive at 
the main result of this paper — reducing the running time to within a polynomial factor 
of o~{G) — requires manipulation at the level of the fast Moebius transform "inside" the al- 
gorithm, which can be found in £4.1 The smooth time-space tradeoff, Theorem [T^c), is 



obtained by a new "split transform" technique (Appendix [B| 

Our approach highlights the algorithmic significance of the Fortuin-Kasteleyn identity, 
and suggests a more general technique: to compute a polynomial, it may be advisable to look 
at its evaluations at integral (or otherwise special) points, with the objective of obtaining 
new combinatorial or algebraic interpretations that then enable faster reconstruction of the 
entire polynomial. (For example, the multiplication of polynomials via the fast Fourier 
transform can be seen as an instantiation of this technique.) 

We also give another vertex-exponential time algorithm that does not rely on interpola- 
tion (£ |4.2[ ). It is based on a new recurrence formula that alternates between partitioning an 
induced subgraph into components and a subtraction step to solve the connected case. The 
recurrence can be solved using fast subset convolution jB] over a multivariate polynomial 
ring. However, an exponential space requirement seems inherent to that algorithm^] Ap- 
pendix [D] briefly reports on our experiences with implementing and running this algorithm; 
it outperforms deletion-contraction in the worst case when n > 13. 

1.3. Conventions. For standard graph-theoretic terminology we refer to West [46J. All 
graphs we consider are undirected and may contain multiple edges and loops. For a graph 
G, we write n = n{G) for the number of vertices, m = m{G) for the number of edges, 
V = V(G) for the vertex set, E = E{G) for the edge set, c = c{G) for the number of 
connected components, t(G) for the number of spanning trees, and cr(G) the number of 
connected sets, i.e., the number of vertex subsets that induce a connected graph. 

To simplify running time bounds, we assume m = and remark that this assumption 
is implicit already in Theorem [TJ (Without this assumption, all the time bounds require an 
additional multiplicative term m 0<yl \) For a set of vertices U C V(G), we write G[U] for 
the subgraph induced by U in G. A subgraph H of G is spanning if V(H) = V(G). For a 
proposition P, we use Iverson's bracket notation [P] to mean 1 if P is true and otherwise. 

2. Prior work: Algorithms for the Tutte Polynomial 

The direct evaluation of Tq{x, y) based on Q takes 2 m n°^ steps and polynomial space, 
but many other expansions have been studied in the literature. 

2.1. Spanning Tree Expansion. If we expand and collect terms in Q we arrive at 
(2) T G (x,y) = Y J ti j x t y J . 

In fact, this is Tutte's original definition. The coefficients of this expansion are well- 
studied: assuming that G is connected, tij is the number of spanning trees of G having 
"internal activity" i and "external activity" j. What these concepts mean need not occupy 
us here (for example, see [U §13]), for our purposes it is sufficient to know that they can 
be efficiently computed for a given spanning tree. Thus @ can be evaluated directly by 
iterating over all spanning trees of G, which can be accomplished with polynomial delay 
|27j . The resulting running time is within a polynomial factor of t(G). 



*"A previous version of this manuscript followed this route, establishing Theorem 1(a) 
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Some of the coefficients tij have an alternative combinatorial interpretation, and some 
can be computed faster than others. For example, ioo = holds if m > 0, and toi = *io 
if m > 1. The latter value, the chromatic invariant 6(G), can be computed from the 
chromatic polynomial, and thus can be found in time 2 n n°^ |8j. 

The computational complexity of computing individual coefficients iy has also been in- 
vestigated. In particular, polynomial-time algorithms exist for t n -\-kj for constant k and 
all j = 0, 1, . . . , m — n + 1. In general, the task of computing t^ is #P-complete PJ. 

2.2. Deletion— Contraction. The classical algorithm for computing T G is the following 
deletion- contraction algorithm. It is based on two graph transformations involving an edge 
e. The graph G\e is obtained from G by deleting e. The graph G/e is obtained from G by 
contracting e, that is, by identifying the endvertices of e and then deleting e. 
With these operations, one can establish the recurrence formula 

1 if G has no edges; 

y T G\e(x,y) if e is a loop; 

xT G / e (x,y) if e is a bridge; 

T G \e(x, y) + T G/e (x, y) otherwise. 

The deletion-contraction algorithm defined by a direct evaluation of ^ leads to a running 
time that scales as the Fibonacci sequence, ((1 + y/b) /2) n+m = O(1.6180 n+m ) [47 . Sekine, 
Imai, and Tani [38j observed that the corresponding computation tree has one leaf for every 
spanning tree of G, so ^ is yet another way to evaluate Tq in time within a polynomial 
factor of r(G). In practice one can speed up the computation by identifying isomorphic 
graphs and using dynamic programming to avoid redundant recomputation \22\ |2~4] [38] . 

The deletion-contraction algorithm is known to compute many different graph parame- 
ters. For example, the number of spanning trees admits an analogous recursion, as does the 
number of acyclic orientations, the number of colourings, the dimension of the bicycle space, 
and so forth |20[ §15.6-8]. This is no surprise: all these graph parameters are evaluations of 
the Tutte polynomial at certain points. But not only is every specialisation of Tq express- 
ible by deletion-contraction, the converse holds as well: every graph parameter that can be 
expressed as a deletion-contraction recursion turns out to be a valuation of Tq, according 
to the celebrated Recipe Theorem of Oxley and Welsh [36] (cf. |10| Theorem X.2]). 

Besides deletion-contraction, many other expansions are known (in particular for restric- 
tions of the Tutte polynomial; see [4]), even a convolution over the set of edges [33 , but 
none leads to vertex-exponential time. 



(3) T G {x,y) = { 



2.3. Regions of the Tutte plane. The question at which points (x,y) the Tutte poly- 
nomial can be computed exactly and efficiently was completely settled in the framework of 
computational complexity in the seminal paper of Jeager, Vertigan, and Welsh [26]: They 
presented a complete classification of points and curves where the problem is polynomial- 
time computable, and where it is #P-complete. This result shows us where we probably 
need to resign ourselves to a superpolynomial-time algorithm. 

For most of the #P-hard points, the algorithms from {2.1 and {2.2 were best known. 
However, for certain regions of the Tutte plane, algorithms running in time exp(0(n)) 
have been known before. We attempt to summarise these algorithms here, including the 
polynomial-time cases; see Figure [T] 
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Trivial hyperbola: On the hyperbola (x — l)(y — 1) = 1 the terms of ([I]) involving 
c{F) cancel, so Tq(x, y) = (x — l) n ~ c y m , which can be evaluated in polynomial time. 

Ising model: On the hyperbola H2 = (x — l)(y — 1) = 2, the Tutte polynomial gives 
the partition function of the Ising model, a sum of easily computable weights over 
the 2 n configurations of n two-state spins. This can be trivially computed in time 
2« n O(i) an( j p ly nom i a l space. By dividing the n spins into three groups of about 
equal size and using fast matrix multiplication, one can compute the sum in time 
2nui/3 n O{i) _ 0(1.732 n ) and exponential space, where u is the exponent of matrix 
multiplication; this is yet a new application of Williams's trick [5j [32j 148] . 

Potts model: More generally, for any integer q > 2, the Tutte polynomial on the 
hyperbola H q = (x — l)(y — 1) = q gives the partition function of the q-state 
Potts model |37| . This is a sum over the configurations of n spins each having q 
possible states. It can be computed trivially in time oTn ^ and, via fast matrix 
multiplication, in time q^/^n ^. We will show in ^3jthat, in fact, time 2 n n°^ 
suffices, which result will be an essential building block in our main construction. 

Reliability polynomial: The reliability polynomial Rc(p), which is the probability 
that no component of G is disconnected after independently removing each edge 
with probability 1 — p, satisfies Rg(p) = p m ~ n+c (l — p) n ~ c Tc(l, 1/p) and can be 
evaluated in time 3 n n°^ and exponential space 

Number of spanning trees: For connected G, Tc(l,l) equals the number r(G) of 
spanning trees, and is computable in polynomial time as the determinant of a maxi- 
mal principal submatrix of the Laplacian of G, a result known as Kirchhoff 's Matrix- 
Tree Theorem. 

Number of spanning forests: The number of spanning forests, Tg(2,1), is com- 
putable in time 2 n n°^ by first using the Matrix- Tree Theorem for each induced 
subgraph and then assembling the result one component (that is, tree) at a time via 
inclusion-exclusion [8 . (This observation is new to the present work, however.) 

Dimension of the bicycle space: Tq(— 1, —1) computes the dimension of the bicy- 
cle space, in polynomial time by Gaussian elimination. 

Number of nowhere-zero 2-fiows: 7g(0, —1) = 1 if G is Eulerian (in other words, 
it "admits a nowhere-zero 2-flow"), and Tg(0, — 1) = otherwise. Thus Tc(0, — 1) 
is computable in polynomial time. 

Chromatic polynomial: The chromatic polynomial Pc(t): which counts the number 
of proper t-colourings of the vertices of G, satisfies Pc(t) = (—l) n ~ c t c TG(l — i, 0) 
and can be computed in time 2 n n°^ [8j. Vertex-exponential time algorithms were 
known at least since Lawler [33], and a vertex-exponential, polynomial-space algo- 
rithm was found only recently |5j. Other approaches to the chromatic polynomial 
are surveyed by Anthony |3_. At t = 2 (equivalently, x = — 1) this is polynomial- 
time computable by breadth-first search (every connected component of a bipartite 
graph has exactly two proper 2-colourings) . The cases t = 3,4 are well-studied 
benchmarks for exact counting algorithms, the current best bounds are 0(1.6262 n ) 
and 0(1.9464 n ) [E]. The case x = is trivial. 

To the best knowledge of the authors, no algorithms with running time exp(0(n)) have 
been known for other real points. If we allow x and y to be complex, there are four more 
points (x, y) at which Tq can be evaluated in polynomial time |26j . 
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2.4. Restricted graph classes. Explicit formulas for Tutte polynomial have been derived 
for many elementary families of graphs, such as T(C n ; x, y) = y + x + x 2 + • • • + x n_1 for 
the n-cycle graph C n . We will not give an overview of these formulas here (see jH §13]); 
most of them are applications of deletion-contraction. 

For well-known graph classes, the authors know the following results achieving exp(0(n)) 
running time or better: 

Planar graphs: If G is planar, then the Tutte polynomial can be computed in time 
exp(0(y / ra)j [38 . This works more generally, with a slight overhead: in classes of 
graphs with separators of size n a , the Tutte polynomial can be computed in time 
exp(0(n Q logn)). 

Bounded tree- width and branch-width: For k a fixed integer, if G has tree-width 
k then Tq can be computed in polynomial time [TJ [35] . This can be generalised to 
branch- width [23 . 

Bounded clique-width and cographs: For k a fixed integer, if G has clique-width 
k then Tq can be computed in time exp {0(n l - l l (k+ ^)) pSJ. A special case of this 
is the class of cographs (graphs without an induced path of 4 vertices) , where the 
bound becomes exp (0(n 2 / 3 )). 

Bounded-degree graphs: If A is the maximum degree of a vertex, the deletion- 
contraction algorithm and 2m < nA yield the vertex-exponential running time 
bound O(l.6180( 1+A / 2 ) n ) directly from the recurrence. Gebauer and Okamoto im- 
prove this to x\ n0( - l \ where xa = 2(1 - A2~ A ) 1 / (A+1 ) (for example, X3 = 2.5149, 
X4 = 3.7764, and X5 = 5.4989). For /c-regular graphs with k > 3 a constant 
independent of n, the number of spanning trees (and hence, within a polyno- 
mial factor, the running time of the deletion-contraction algorithms) is bounded 
by t(G) = 0(i^n _1 log n), where u/~ = (k — l) k ^ 1 /(k 2 — 2k) k / 2 ~ 1 (for example, 
u 3 = 2.3094, z/ 4 = 3.375, and u 5 = 4.4066), and this bound is tight [H]. 

Interval graphs: If G is an interval graph, then Tq can be computed in time 
O(1.9706 m ), which is not exp(0(ra)) in general, but still faster than by deletion- 
contraction |17| . 

What we cannot survey here is the extensive literature that studies algorithms that 
simultaneously specialise Tq and restrict the graph classes, often with the goal of developing 
a polynomial-time algorithm. A famous example is that for Pfaffian orientable graphs, which 
includes the class of planar graphs, the Tutte polynomial is polynomial-time computable 
on the hyperbola H2 [29 . Within computer science, the most studied specialisation of this 
type is most likely graph colouring for restricted graph classes. 

2.5. Computional complexity. The study of the computational complexity of the Tutte 
polynomial begins with Valiant's theory of #P-completeness [42 j and the exact complexity 
results of Jaeger, Vertigan, and Welsh [26]. The study of the approximability of the values 
of Tq has been a very fruitful research direction, an overview of which is again outside the 
scope of this paper. In this regard we refer to Welsh's monograph [13] and to the recent 
paper of Goldberg and Jerrum |21] for a survey of newer developments. 

For our purposes, the most relevant hardness results have been established under the 
Exponential Time Hypothesis [25] (ETH). First, deciding whether a given graph can be 
3-coloured requires exp(0(n)) time under ETH, and since 3-colourability can be decided 
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by computing T G (— 2,0) we see that evaluating the Tutte polynomial requires vertex- 
exponential time under ETH. Thus, it would be surprising if our results could be significantly 
improved, for example to something like exp (0(n/ log n)). 

Second, it is by no means clear that the entire Tutte plane should admit such algorithms. 
Many specialisations of the Tutte polynomial can be understood as constraint satisfaction 
problems. For example, graph colouring is an instance of (g, 2)-CSP, the class of constraint 
satisfaction problems with pairwise constraints over (/-state variables. Similarly, the par- 
tition function for the Potts model can be seen as a weighted counting CSP [19]. Very 
recently, Traxler [4(D has shown that already the decision version of (g, 2)-CSP requires 
time exp (r2(nlogg)) under ETH, even for some very innocent-looking restrictions, and 
even for bounded degree graphs. Thus in general, these CSPs are not vertex-exponential 
under ETH. 

3. The multivariate Tutte polynomial via the (/-state Potts model 

Let R be a multivariate polynomial ring over a field and let G be an undirected graph 
with vertex set V = {1,2, . . . , n} and edge set E, m = n 0<yl \ We allow G to have parallel 
edges and loops. Associate with each e 6 E a ring element r e G R. The multivariate Tutte 
polynomial [39] of G is the polynomial 

(4) Z G (q,r)= 5> c(F) n^> 

FCE eG-F 

where q is an indeterminate and c(F) denotes the number of connected components in the 
graph with vertex set V and edge set F. The product over an empty set always evaluates 
to 1. 

The classical Tutte polynomial T G {x, y) can be recovered as a bivariate evaluation of the 
multivariate polynomial Z G (q, r) via 

(5) T G (x,y) = (x - l)~< E \y - l)-^Z G ((x - l)(y - l),y - l) . 

3.1. The Fortuin Kasteleyn identity. At points q = 1,2,... the multivariate Tutte 
polynomial Z G (q,r) can be represented as an evaluation of the partition function of the 
g-state Potts model [161 EH] • 

For a mapping s : V — > {1, 2, . . . , q} and an edge e £ E with endvertices x and y, define 
8% = 1 if s(x) = s(y) and 8% = if s(x) / s(y). The partition function of the g-state Potts 
model on G is defined by 

(6) Z^(q,r)= Y, n( 1+r ^)- 

s:V->{l,2,...,g} e£E 

Theorem 3 (Fortuin and Kasteleyn). For all q = 1, 2, . . . it holds that 

(7) Z G (q,r) = Zg**(q,r). 

3.2. The multivariate Tutte polynomial via the g-state Potts model. By virtue of 
the Fortuin-Kasteleyn identity ([7]), to compute Z G (q,r) it suffices to evaluate 

Zg otts (l,r), Zg otts (2,r), Zg otts (n + 1, r) 

and then recover Z G (q, r) via Lagrangian interpolation. For the interpolation to succeed, 
it is necessary to assume that the coefficient field of R has a large enough characteristic so 
that 1,2, ... ,n have multiplicative inverses. 
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At first sight the evaluation of ^ for a positive integer q appears to require q n n°^ ring 
operations. Fortunately, one can do better. To this end, let us express ZQ° tts (q, r) in a more 
convenient form. For icy, denote by G[X] the subgraph of G induced by A, and let 

(8) f{x)= n 

eeE(G[X]) 

For q = 1, 2, . . ., we have 

(9) Z^(q,r) = Y, Wi)f(U 2 )---f(U g ), 

(Ul,U2,...,Ug) 

where the sum is over all g-tuples (U\, U2, ■ ■ ■ , U q ) with U\, U2, ■ ■ ■ , U q C V such that 
U q i=1 Ui = V and Uj n t7 fc + for all 1 < j < k < q. 

We now proceed to develop algorithms for evaluating the Potts partition function in the 
form @. 

3.3. The baseline algorithm. Let / : 2 V — > i? be a function that associates a ring element 
/(A) G i? with each subset X (IV. 

The zeta transform f( : 2 V -> R is defined for all F C V by /C(F) = J2xcy f( x )- The 
Moeferas transform fp : 2 y -»■ is defined for all X C V by /^(X) = EycxC" 1 ) Wl 

It is a basic fact that the zeta and Moebius transforms are inverses of each other, that 
is, /£/z = = / for all /. Furthermore, it is known [6 that 

(10) ((/C)V)W= E 

(E/i,f/ 2j ...,f/ s ) 

where the sum is over all g-tuples (Z7i, t/2, • • • , C^g) with C7i, C/2, ...,U 9 CF and Uj =1 Uj = V. 
In particular, ((/C) 9 ^)(^) can t> e computed directly in 3 n n°W ring operations by storing 
n O(i) r j n g e l emen ts. Using the fast zeta and Moebius transforms, ((/C) 9 /^) (V) can be 
computed in 2 n n°W ring operations by storing 2 n n°^ ring elements [6]. 

To use this to evaluate Q, adjoin a new indeterminate z into -R to obtain the polynomial 
ring R[z]. Replace / with f z : 2 V -» defined for all X C V by / Z (X) = /(X)zl x L Now 
evaluate the z-polynomial ((f z () q fjMV) and look at the coefficient of the monomial 
which by virtue of (10) is equal to (J9|. 



This baseline algorithm together with (|5]), ([7]), and Lagrangian interpolation establishes 
that the Tutte polynomial Ta(x,y) can be computed (a) in time and space 2 n n°^; and 
(b) in time 3 n n°^ and space n°^\ This proves Theoremjlm) . A more careful analysis of 
((fC) q /j)(V) enables the time-space tradeoff in Theorem j^c). [[ See Appendix [b| ]] 

4. Improvements and variations 

4.1. An algorithm over connected sets. It is useful to think of A C V in what follows as 
the current subset under consideration. We start with a lemma that partitions the subsets 
of X based on the maximum common suffix. To this end, let Y=% A be a shorthand for 
Y n {i + l,i + 2, . . . ,n} = X n {i + l,i + 2, . . . ,n}. 

Lemma 4 (Suffix partition). Let Y C A C {1, 2, . . . , n}. Then, either Y = X or there 
exists a unique i G A such that Y=i- \ X \ {i}. 

Proof. Either Y = X or i = max A \ Y. □ 
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The intermediate values computed by the algorithm are now defined as follows. 
Definition 5. Let X C V, q = 1, 2, . . . , n + 1, and i = 0, 1, . . . , n. Let 

F(X,q,i)= 11^)' 
(Ui,U 2 ,...,U q ) j=l 

where the sum is over all g-tuples (Ui, U2, • • • , U q ) such that both J/j, U2, ■ ■ ■ ,U q C. X and 

Note that F(V,q,0) = ((fO q ^)(V). Thus, it suffices to compute F(V,q,0). 

We are now ready to describe the algorithm that computes the intermediate values 
F(X, q, i) in Definition |5j The algorithm considers one set X C V at a time, starting 
with the empty set X = and proceeding upwards in the subset lattice. It is required that 
the maximal proper subsets of X have been considered before X itself is considered; for 
example, we can consider the subsets of V in increasing lexicographic order. The comments 
delimited by "[[" and "]]" justify the computations in the algorithm. 

Algorithm U. (Up-step.) Computes the values F(X,q,i) associated with X using the 
values associated with X \ {i} for all i G X. 

Input: A subset icy and the value F(X\{i}, q, for each i £ X and q = 1,2, . . . , n+1. 
Output: The value F(X, q, i) for each q = 1, 2, . . . , n + 1 and i = 0, 1, . . . , n. 
Ul: For each q = 1, 2, 3, . . . , n + 1, set 

q, n) = (f{X) + ^F(X\ {i}, 1, * - L ^ 

[[ By the suffix partition lemma, J2ycx /00 = Siex \ W> 1? * — !)• Adding 
f(X) and taking powers, we obtain F(X,q,n). ]| 

U2: For each g = 1, 2, 3, . . . , n + 1 and i = n, n — 1, . . . , 1, set 

F(X, g, i - 1) = F(X, q, i) - [i G \ {t}, g, i - 1) . 

[ There are two cases to consider to justify correctness. First, assume that i £ X. 
Consider an arbitrary g-tuple (U\, U2, ■ ■ ■ , U q ) with Ui, f7 2 , . . . , U q C X. Let Y = 
U Q j =l Uj. Clearly, Y C X. Because i X and Y C X, we have K=i_i X if and 
only if y=j X. Thus, F(X,q,i — 1) = F(X,q,i). Second, assume that i G X. In 
this case we have y=j X if and only if either Y=j_i A" or Y=i-\ X\{i} (the former 
case occurs if i G Y, the latter if i ^ y). In the latter case, Y C A \ {i} and hence 
U x ,U 2 ,... l U q ^X\ {%}. Thus, g, i - 1) = F(A, q, i) - F(X \ {i}, g, i - 1). ]] 

Assume that / satisfies the following property: for all X C V it holds that 

(11) f(X) = f(X 1 )f(X 2 ).-.f(X 8 ) 

where G[A"i], G[X2], . . . , G[A" S ] are the connected components of G[X]. For convenience 



we also assume that /(0) = 1. Note that the factorisation (11) is well-defined because of 



commutativity of R. Also note that <\8k satisfies (11) 



Lemma 6. Let Gr[Xi], GfXa], . . . , Gf-Xg] be the connected components of G[X\. Then, 

s 

(12) F(X,q,i) = llF(X k ,q,i). 



k=l 
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The recursion (12) now enables the following top-down evaluation strategy for the inter- 
mediate values in Definition |5j Consider a nonempty X C V . If G[X] is not connected, 
recursively solve the intermediate values of each of the vertex sets X\ , X2 , • • • , X s of the 
connected components GLXi], GLY2], . . . , G[A" S ] of G[X], and assemble the solution using 
(12). Otherwise; that is, if G[X] is connected, recursively solve the intermediate values of 
each set X \ {i}, i S X, and assemble the solution using Algorithm U. Call this evaluation 
strategy Algorithm C. 

Algorithm C together with (|5]), ([7]), and Lagrangian interpolation establishes that the 
Tutte polynomial T G (x,y) can be computed in time and space a{G)n°^ 1 \ This proves 
Theorem [TJa). 

4.2. An alternative recursion. We derive an alternative recursion for Z G (q,r) based on 
induced subgraphs and fast subset convolution. Let R be a commutative ring. Associate a 
ring element r e £ R with each e 6 E. For k = 1, 2, . . . , n, let 



S G (k,r)= yi n 



r,, 



FCE eSF 
c(F)=k 

and observe that Z G (q,r) = ^2k=i1 k ^G(k,r). Thus, to determine Z G (q,r), it suffices to 
compute S G (k, r) for all k = 1, 2, . . . , n. 

To this end, the values Sc(k,r) can be computed using the following recursion over 
induced subgraphs of G. Let W C V and consider the subgraph G[W] induced by W in G. 
Suppose that SQp](k, r) has been computed for all 7^ U C^W and k = 1, 2, . . . , \ U\. 

To compute S G [ W ](k, r) for k = 2, 3, ... , \W\, observe that a disconnected subgraph of 
G[W] partitions into connected components. Thus, for k > 2 we have 

(13) S G iw](k,r) = - Y Sciui^^^S^w^ik -l,r) . 

For the connected case, that is, for k = 1, it suffices to observe that we can subtract the 
disconnected subgraphs from the set of all subgraphs to obtain the connected graphs; put 
otherwise, 

(14) S G[w] (l,r)= J] (l + r e )-J2s G [w](k,r). 

e&E(G[W]) k>2 



The recursion defined by (13) and (14) can now be evaluated for \W\ = l,2,...,n in 
total 2 n n°( 1 ) ring operations using fast subset convolution [B]. As a technical observation 
we remark that (13) assumes that k has a multiplicative inverse in R; this assumption 
can be removed, but we omit the details from this extended abstract. We also note that 
analogues of Algorithms U and C running in a(G)n°^ ring operations can be developed 
in this context; we describe an implementation of this in Appendix [D] However, it is not 
immediate whether a polynomial-space algorithm for the Tutte polynomial can be developed 
based on (fl~3l) and (fl4|). 
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Appendix 

Appendix A. Proofs 

A.l. Proof of Theorem |3j This proof of the Fortuin-Kasteleyn identity ([7]) is well known 
(e.g. |39j ) and is here included only for convenience of verification. 

Proof. Expanding the product over E and changing the order of summation, 

z p G otts ( q ,r)= e II( 1 +^e) = E E II^- 

s:V->{l,2,...,g} e£E FCE s:V—>{l,2,...,q} eCF 

The right-hand side product evaluates to zero unless s is constant on each connected com- 
ponent of the graph with vertex set V and edge set F. Because there are q choices for the 
value of s on each connected component, 

e e n^=E? c(F) n^=^(^). 

FCE s:V->{l,2,...,q}e<=F FCE e£F 

□ 

A. 2. Proof of Lemma [6[ It is convenient to start with a preliminary lemma. 

Lemma 7. Let G[Xi],G[X2], ■ ■ ■ ,G[X S ] be the connected components of G[X] and let U C 
X. Then, 

f(U) = f(u n x x )f{u ni 2 )- f(u n x.) . 



Proof Let G[Ui],G[U2], ■ ■ ■ ,G[Ut] be the connected components of G[U]. Then, by (11), 

f(U) = f(U 1 )f(U 2 )---f(U t ). 

Because PCI holds, for every Ui there is a unique h(i) E {1, 2, . . . , s} such that Ui C X h ^y 
Moreover, since {Ui, U2, ■ ■ ■ , Ut} is a partition of U, we have that {Ui : % G is a 

partition of \JV\Xj for all j = 1, 2, . . . , s. Thus, by fy]) we have f(UC\Xj) = Uieh-^(j) f( U i) 



for all j = 1, 2, . . . , s. In particular, by commutativity of R, 

t s s 

f(u) = Uf(u t ) = H J] f(u l ) = l\f(unx J 

i=l J^lie/i-iQ') j=l 



□ 



We now proceed with the proof of Lemma [6j 
Proof. Consider an arbitrary g-tuple (Ui, U2, ■ ■ ■ , U g ) with U±, U2, ■ ■ ■ , U q CI and yj q , =l Uj=i X. 



Because {X±, X2, ■ ■ ■ ,X S } is a partition of A, we have U q - =1 Uj=iX if and only if X^ n 



r- 

U q =1 Uj=i Xk n X holds for all k = 1, 2, . . . , s. Put otherwise, we have \J q =1 Uj=i X if and 
only if uJ =1 (Xfc n Uj)=i X^ holds for all k = 1, 2, . . . , s. Using Lemma [7] for each Uj in turn, 
we have, by commutativity of R, the unique factorisation into pairwise intersections 



/(*7i)/(c/ 2 ) . . . nu q ) = n n Wj n = n n /(^ n ^) • 

jr = lfc=l fc=lj = l 

The claim follows because (C/i, U2, ■ ■ ■ , C/q) was arbitrary. □ 
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Appendix B. A time-space tradeoff via split transforms 

This appendix outlines a "split transform" algorithm that enables a time-space tradeoff 
in evaluating ((/C) 9 /^) (Y) f° r a given function / : 2 V — ► R and q = 1, 2, . . . , n + 1. 

Split the ground set V = {1,2, ... , n} into two parts, V\ C y and V2 C V, such that 
y = Vi U y 2 and Vi n V2 = 0. Let m = |Vi| and n 2 = \V 2 \. For a subset IC^,we use 
subscripts to indicate the parts of the subset in V\ and V 2 ; that is, we let X\ = X C\V\ and 
X 2 = X PI V 2 . It is also convenient to split the function notation accordingly, that is, we 
write f(Xi,X 2 ) for f(X\ UI2) = f(X). In the context of zeta and Moebius transforms, we 
use X for a subset in the "spatial" (original) domain and Y for a subset in the "frequency" 
(transformed) domain. 

An elementary observation is now that both the zeta and Moebius transforms split, that 

is, 

/coo=£/po = E E f( x ^ x ^) = E ^p^H^ac^) 

XCY X 1 CY 1 X 2 CY 2 XiQYi 

and 

Mi)=^(-i)^/(y)= £ ( _i)PWi £ (-1)^1/(7,7,) 

YCX ATiCYx X2CY2 

= Y, (-l) lXAYl] f»2(Y u X 2 ) = fwi(X 1 ,X 2 ). 

XiQYi 

Also note that /( = /C2C1 = /C1C2 and f/j, = fn 2 /J,i = fmn 2 . 

To arrive at the split transform algorithm for computing ((f() q /J,) (V), split the outer 
Moebius transform and the inner zeta transform to get 

((/ov)(to= E (-i) imYi1 E (-i) imY2l (/CiC2(7i,7 2 ))^ . 

Now let Yi be fixed and consider the inner sum. To evaluate the inner sum for a fixed Y±, 
it suffices to have f(i( 2 (Yi,Y 2 ) available for each Y 2 C V 2 . By definition, 

fCMY 1 ,Y 2 ) = Y, Ki(Yi,X 2 ). 
x 2 cy 2 

Observe that if we have fCi(Yi, X 2 ) stored for each X 2 C V 2 , then we can evaluate 
f(i(2(Yi,Y 2 ) for each Y 2 C V 2 simultaneously using the fast zeta transform. This takes 
in total at most 2 n2 n 2 ring operations and requires one to store at most 2 n2 n 2 ring ele- 
ments. 

For fixed Y\ and X 2 , we can evaluate and store 

f( 1 (Y 1 ,X 2 )= Y f(Xi,X 2 ) 

XiCYi 

by plain summation in at most 2'^' ring operations. Thus, for fixed Y\, we can evaluate 
fC,i{Yi,X 2 ) for each X 2 C V 2 in total at most 2l Yl l2™ 2 ring operations. 

Considering each Y\ C V\ in turn, we can thus evaluate ((/C) 9 ^) (V) by storing at most 
2 n2 n 2 ring elements and executing at most 

n O(i) £ (2" 2 n 2 +2l yi l2 n2 ) = n°( 1 )(3 ril +2 ni n 2 )2 n2 

ncvi 
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ring operations. This completes the description and analysis of the split transform algo- 
rithm. 

The split transform algorithm together with ^ , ([7| , and Lagrangian interpolation proves 
Theorem [ljc). 

Appendix C. The cover polynomial 

Let D be a digraph with vertex set V = {1,2,. . . , n}. Note that D may have parallel 
edges and loops. We assume that the number of edges is n°^ l \ Denote by cr>(i,j) the 
number of ways of disjointly covering all the vertices of D with % directed paths and j 
directed cycles. The cover polynomial is defined as 

C D (x,y) = ^2c D (i,j)x l -y J , 

where x- = x(x — 1) • • • (x — i + 1) and x- = 1. It is known that Cd(x, y) is #P-complete 
to evaluate except at a handful of points (x, y) [9]. 

In analogy to Theorem [TJ we can show that Cjj can be computed in vertex-exponential 
time: 

Theorem 8. The cover polynomial of an n-vertex directed graph can be computed 

(a) in time and space 2 n n°^ ; and 

(b) in time 3 n n°^ and polynomial space. 

The proof involves several inclusion-exclusion-based arguments with different purposes 
and in a nested fashion, so we first give a high-level overview of the concepts involved. One 
readily observes that the cover polynomial can be expressed as a sum over partitionings 
of the vertex set, each vertex subset appropriately weighted, so the inclusion-exclusion 
technique [8] applies. Computing the weights for all possible vertex subsets is again a hard 
problem, but the fast Moebius inversion algorithm [7] can be used to compute the necessary 
values beforehand. This leads to an exponential-space algorithm. Finally, to use inclusion- 
exclusion to reduce the space to polynomial [28, I5T] , we apply the mentioned transforms in 
a nested manner and switch the order of certain involved summations. 

We turn to the details of the proof. For Xcy, denote by p{X) the number of spanning 
directed paths in -D[X], and denote by c{X) the number of spanning directed cycles in 
D[X\. Define p(0) = c(0) = 0. Note that for all x G V we have p{{x}) = 1 and that c({x}) 
is the number of loops incident with x. 

By definition, 

c D (hj) = ~ p(X 1 )p(X 2 )---p(X i )c(Y 1 )c(Y 2 )---c(Y j ), 

% ' 3 ' X 1 ,X 2 ,...,X l ,Yi,Y2,-.,Yj 

where we sum over all (i+j)-tuples (Xi, X%, . . . , JQ, Y±,Y2, ■ ■ . , Yj) such that {X%, X 2 , ■ ■ ■ , Xi, 
Yi, Y 2 , . . . , Yj} is a partition of V. 

We next derive an alternative expression using the principle of inclusion and exclusion. 
To this end, it is convenient to define for every f/C^ the polynomials 



P(U;z) = ^ p{X)z\ x \ and C{U;z) = ^ c{X) Z \ x \ 
xcu xcu 
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in an indeterminate z; if viewed as set functions, P(U ; z) and C(U; z) are zeta transforms 
of the set functions p(X)z^ x \ and c{X)z )[X \ respectively. We can now write 

ucv 

It remains to show how to compute the p{X) and c{X) for all X C V. For 5" C V let 
w(S,s,t,£) denote the number of directed walks of length I from vertex s to vertex t in 
D[S]; define w(S, s, t, I) = if s S or t 5. By inclusion-exclusion, again, 

Ks<KnSCX 

Similarly, 

c{X)= ^2(-l)^ x \ s ^w{S,s,s, \X\), where s = minS. 

sex 

Observing that w(S,s,t,£) can be computed in time n°^ l \ we have that C[>(i,j) can be 
computed in space and time A n n°^ 1 \ 

To get an algorithm running in 3 n n ^ time and space, observe that 

P(U;z)=Y,P(U,S;z) 

scu 

where 

P(U,S;z)= £ ^( |l/ \ 5| )(-l)^l+ fc ti,(5, a ,t,|5| + fc-l) 

l<s<t<n fc=0 ^ ' 

and 

C(C/;z)= ^C(C/,5;z) 
scu 

where 

|CAS| /ir/\ <?i\ 

C{U,S;z)= r ^ l ){-l) k z lsl+k w{S,s,s,\S\ + k), where s = minS. 
k=o ^ ' 
This establishes part (b) of the theorem. 

For part (a), we show how to evaluate cd in time and space 2 n n°( 1 ). Namely, p and 
c can be computed in time and space 2 n n°^ via fast Moebius inversion. Given p and c, the 
polynomials P and C can be computed in time and space 2 n n°^ via fast zeta transform. 
And finally, given P and C, the inclusion-exclusion expression of cu(i,j) can be evaluated 
in time 2 n n° ( - 1 \ 

Appendix D. Tutte polynomials of concrete graphs 



D.l. Algorithm implementation. Our implementation of the algorithm described in {4.2 



uses a number of extra techniques to reduce the polynomial factors in the time and memory 
requirements. In what follows we assume that G is a connected graph. 

(1) The coefficients Uj of the Tutte polynomial are computed modulo a small integer 
p; the computation is repeated for sufficiently many different (pairwise coprime) 
p to enable recovery of the coefficients via the Chinese Remainder Theorem. The 
number of different p required is determined based on the available word length and 
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using t{G) (computed via the Matrix- Tree Theorem) as an upper bound for the 
coefficients. 

(2) To save a factor of m in memory, instead of direct computation with bivariate 
polynomials, we compute with univariate evaluations of the polynomials at z = 
0,1, ... ,m, and finally recover only the necessary bivariate polynomials from the 
evaluations via Lagrange interpolation. 

(3) To save a further factor of n 2 in memory, we execute the analogue of Algorithm U 
for subsets X in a specific order, namely in the lexicographic order. This enables 
efficient "in-place" computation of the polynomials F(X, k, i) so that, for each X, 
the polynomials F(X, k, i) need to be stored only for one value of i at the time. 
Furthermore, we never need all F{X, k, i) for k = 2, 3, . . . , n explicitly, only a linear 
combination of them, so we count with this instead; however, we omit the details in 
this abstract. 

The source code of the algorithm implementation is available by request. The imple- 
mentation uses the GNU Multiple Precision Arithmetic library (http://gmplib.org/) for 
computation with large integers. The computed coefficients Uj are checked for consistency 
by verifying that Uj = t{G) and that J2i,j 2 i+i % = 2 m . 

D.2. Performance. The current algorithm implementation uses roughly 2 n+1 n words of 
memory for an n- vertex graph, which presents a basic obstacle to practical performance. For 
example, the practical limit is at n = 25, assuming 32 GB of main memory and 64-bit words. 
This makes our polynomial space and time-space tradeoff algorithms from Theorem [T](b,c) 
interesting also from a practical perspective. At the time of writing, we have implemented 
the former, but not yet performed large-scale experiments with it. 

In terms of running time, the complete graph K n presents the worst case for n-vertex 
inputs for our algorithm. On a 3.66GHz Intel Xeon CPU with 1MB cache, computing the 
Tutte polynomial of Kn takes less than an hour, K\% takes about three hours, and K22 
takes 96 hours. In comparison, both deletion-contraction and spanning tree enumeration 
cease to be practical well below this; for example, t(K 22 ) = 705429498686404044207947776 
and t(K\q) = 72057594037927936; a survey of how to compute Tq in practice [23] reports 
running times for the complete graph K14 in hours. The fastest current program to compute 
Tutte polynomials [22] is also based on deletion-contraction with isomorphism rejection, but 
uses many other ideas as well. It processes K14 and many sparse graphs with far larger n 
in a few seconds, but also ceases to be practical for some dense graphs with n = 16, see 
Figure [2} 

Two further remarks are in order. First, for (connected) graphs with a small t(G), 
enumeration of spanning trees is faster than our algorithm. Second, graphs with fewer 
edges are faster to solve using our algorithm. For example, a 3-regular graph on 22 vertices 
can be solved in about five hours. 

D.3. Tutte polynomials of some concrete graphs. Even though few readers are likely 
to derive any insight from the fact that the coefficient of x 2 y 2 in the Tutte polynomial of 
Loupekine's Second Snark is 991226, we feel it germane to our paper to actually compute 
some Tutte polynomials. We include tables of the nonzero coefficients Uj in the expan- 
sion ^ for a number of graphs. Among these, the values for the Petersen graph are well 
known §13b] and are included here for verification only. For reference, we present the 
Tutte polynomials of a few other well-known graphs, mostly snarks and cages; however, 
these graphs are fairly sparse and exhibit symmetries that make them amenable to many of 
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Figure 2. Running times for complements of random 4- regular graphs. The 
lines show averages of 5 runs on a 3.66GHz Intel Xeon CPU with 1MB cache. 
The thin line is our algorithm; the thick line is the algorithm of Haggard, 
Pearce, and Royle |22j . 



the previously existing techniques. An entertaining example that tests the liminations of our 
current implementation is from Knuth's Stanford Graph Base [30 , based on the encounters 
between the 23 most important characters in Twain's Huckleberry Finn. This graph has 
23 vertices, 88 edges, and 54540490752786432 spanning trees; the required solution time is 
about 50 hours. 
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122 002 


38785 


562 021 


233 095 


45 170 


11 144 


1895 


485 


40503 


8615 


704 


144 






2 061 


150 











9 456 659 
58 691487 
163 670 254 
274 803 257 
312 764 825 
256 878 631 
157 948 712 
74 347 451 
27 140739 
7 732 248 
1720 587 
297 808 
39723 
4 020 
301 
16 



j=6 
68 488 
844 000 
561 170 
750421 
917711 
331512 
034 296 
684 553 
158 980 
143 990 
723 236 
086 966 
622 606 
794 424 
529 428 
128 290 
575 231 
11796 



103 243 
264 694 
409 580 
429 835 
325 285 
183 994 
79489 
26 556 
6900 
1394 
218 
26 



j=r 

074 885 356 
788 227 460 
234 344 442 
434 157 718 
626 654 435 
245 146 851 
333 180 127 
180 483 326 
188 103 022 
191730 740 
731774 739 
231 222 093 
157 204 332 
360 083 348 
155 871 340 
7184 647 
210 107 
3005 



30 940 
162 102 
385158 
553 143 
538 846 
378230 
198 149 
79122 
24 368 
5818 
1076 
153 
16 



731400 736 
328114 667 
972 944 741 
680 550 286 
162 737 722 
381 624 169 
002 209146 
703 627 690 
441 736 674 
100 553 716 
252 268 769 
327 979 453 
622 431041 
344 595 660 
78 646 816 
3151857 
77 523 
862 



47 560 
231 997 
514 385 
689 925 
627 608 
411029 
200 617 
74 482 
21270 
4 692 
798 
104 
10 



610 835 303 
488 723 635 
945 146 530 
324 841 686 
139 462 721 
358 220 644 
849105 619 
208 530 196 
851 688 207 
871669 682 
621154 805 
075 037 503 
246 508 473 
745 813 853 
38 768 901 
1 352 207 
27 469 
206 



67 267 401 
307 596 180 
640 213 989 
806 372 976 
688 582 006 
422 898 675 
193 255 210 
67 024 385 
17 826 935 
3 648 698 
573 141 
68506 
6137 
402 
18 



3=10 
905 395 
791466 
471 637 
235 758 
779 197 
284 303 
084 387 
760 956 
293 467 
696 018 
852 737 
947 878 
792 945 
557 913 
595 383 
561013 
9 090 
33 



88 747 
382 634 
751 445 
893 029 
719045 
415871 
178 629 
58078 
14 431 
2 747 
399 
43 



842 449 432 
254 447 952 
805 849 377 
971410 258 
082 827 230 
654 487 327 
421 459 750 
504 339 287 
482 593 462 
096198 020 
057 840 816 
803 197 182 
573 570 819 
211208 919 
8 660 607 
224 483 
2 783 



110 413 
451 056 
839 500 
945 126 
720 198 
393 582 
159 375 
48 701 
11328 
2 008 
270 



J =12 
663 746 802 
497 992 889 
833 190 661 
777 109 455 
797 592 595 
445 653 266 
433 931029 
630 168 789 
241057 760 
371282 725 
019 505 090 
225 629 630 
022 552 038 
107 776 845 
3 932 249 
88116 
859 



130 666 
507 874 
899 138 
962 134 
695 912 
360 278 
137 832 
39 652 
8643 
1428 
177 
16 



= 13 
491 664 749~~ 
432 372 429 
156 074 953 
122 445 120 
053117113 
181 979 195 
995 815 031 
159 130 588 
948 300111 
133 964 932 
714 498 469 
453312 517 
112 678 888 
53 497 484 
1 739 850 
33483 
234 



148 114 
549 640 
928 434 
946 802 
651526 
320 142 
115885 
31417 
6 421 
988 
113 



138 481 
378551 
463 156 
335 152 
159 540 
52 529 
12 241 
2 032 
239 



19 



35 321028 
72 405 228 
64 373 528 
32 605 441 
10 370 366 
2149 485 
291 838 
25 513 
1386 
44 



3=14 
503 232 862 
583 945 429 
404 732 381 
956 316 230 
551686197 
255 061978 
591893 630 
553 021990 
813 488 934 
695 612 555 
804 068 059 
665 564 617 
594 641211 
25 811941 
748 818 
12 227 
53 
3=22 
043 916189 
418 825180 
138 345 255 
391009 976 
234 730 431 
965 418 036 
838 274 716 
190 992 977 
034 177 695 
668 089 937 
118 570 285 
44 326 408 
1 276 909 
25927 
231 
3=30 
354 267 
468100 
302 495 
633781 
228 563 
588 138 
357 248 
055 441 
738494 
763 707 
831 602 
9 466 
72 



161710 933 
574 575 879 
928 302 660 
904 076 789 
592 898 453 
276 885 385 
94 920 219 
24 262 330 
4 649 813 
666 682 
70 899 
5 515 
308 
12 



j =15 
451019 
264 822 
980 981 
608 510 
733529 
561 807 
709713 
193 502 
570718 
299 140 
825 376 
962 845 
435 706 
095 111 
312 758 
4 231 



124191 
327 779 
385 979 
267 823 
121698 
38 034 
8 353 
1294 
140 
10 



3=23 
059 706 563 
450 801 827 
736157 458 
161239724 
403 452 623 
479 409 317 
299 039 147 
966 296 360 
635158 743 
542 647 994 
539 500 322 
19181 752 
502 053 
9 067 
56 

j=31 



27 569 839 
54 377 952 
46 281 448 
22 303 581 
6 697 239 
1 297 255 
162 366 
12 846 
616 
16 



915 839 
292 000 
067 004 
464 318 
653 980 
575 747 
310 290 
282 872 
213381 
929 544 
252 202 
2112 
11 



170813 
582 448 
901802 
840 089 
525 705 
233 534 
75855 
18280 
3 283 
437 
42 



j=16 
927 289 507 
562 263 319 
324 461765 
439 442 061 

387 527 193 

388 363 802 
471 078 622 
905 313 399 
163 455 432 
890 395 312 
952 342 270 
055 400 762 
155 158 798 

5 504 691 
126 758 
1375 



109 337116 
278 622 686 
315 714 271 
209 968 589 
91 005 190 
26 964153 
5570757 
804 342 
80349 
5 459 
249 



= 24 
734795" 
027 317 
443 996 
993 354 
493 607 
070 620 
950 696 
907 429 
078 122 
753 151 
848 063 
926 985 
187 522 
2 949 
10 

= 32 



175 183 
574 307 
853 408 
761 337 
454 

192 359 
59 206 
13447 
2 261 
280 
25 



j =17 
508 018 224 
823 880 181 
768 453 265 
884 261 647 
478155 838 
454 489 345 
435 569 975 
888 149 463 
157 582 200 
092 719 250 
286 305 141 
641278 936 
75 679 090 
2 437 347 
50 082 
428 



94 549 658 
232 599 415 
253 533 343 
161514 
66 709 745 
18 711968 
3629 009 
486 565 
44 519 
2 726 
110 



j =25 
362 427 
134 395 
091 830 
952 720 
029 275 
201 815 
765 916 
736 128 
508 330 
196 170 
746 672 
111957 

65933 
883 



174 936 105 
552153 562 
788 339 037 
674 053 687 
384 901878 
154 895 824 
45166410 
9 661 849 
1518926 
174390 
14453 
854 
35 



j=is 

230 968 
397 919 
061 129 
878 178 
667 513 
102 503 
690 866 
230 137 
901479 
514 342 
004 774 
165572 
781 667 
050 538 
19202 
121 



-26 



80 342 944 
190 758 701 
199 919 236 
121905 869 
47 926 985 
12 706 067 
2 307 827 
286 375 
23887 
1309 
46 



707 533 
968 316 
976 320 
974 532 
451 590 
491 372 
423 945 
736913 
371 320 
543395 
804 725 
152 202 
21528 
239 



170 475 040 
518 597 408 
711981 542 
583 766 784 
318 587 168 
122 026 125 
33 691 771 
6780 631 
994 941 
105615 
8011 
430 
16 



3=19 
838 235 
034 174 
551881 
630 227 
246 233 
706 545 
782 446 
931 834 
569 885 
119 324 
736 241 
162 989 
390 015 
440 428 
7108 
31 



3=27 

67 106 067181304 
153717 226 791019 
154 800 992 492 735 
90 270 537 601981 
33 737 964 980 454 
8 438 244 908 275 
1431 651832106 
163 814154 069 
12 391822 214 
603 613124 
18 772 999 
398 423 
6 401 
56 



162 410 
476 549 
629 449 
495 042 
258 204 
94 096 
24 580 
4 647 
635 
62 



j =20 
989 874 779 
124 971777 
797 206 205 
677 879 502 
924 355 942 
567 995 260 
817 533 069 
801211901 
198 727 558 
160 275 509 
301315057 
209 314 543 
7 263 816 
178942 
2 503 



j =28 

55104 744 824 358 
121719136 468 979 
117 698 380 405 855 
65 567 617176 755 
23 261782 819 227 
5 477 591996 969 
865 610 892 679 
90 962 163 938 
6 204115029 
266 217 054 
7105935 
127 018 
1684 
10 



151482 
428 948 
545 264 
411 367 
205 018 
71045 
17 540 
3110 
394 
35 



3=21 
203 761928 
242 539 842 
254 521678 
624 748 301 
507 541 185 
733 735 891 
994 280 576 
824 903 986 
989 616 510 
513 089 086 
232 878 866 
98 214 601 
3 106 323 

69 836 



44491423252 883 
94 708 476115 968 
87 858 171 960160 
46 701038 037 329 
15 701786 388 350 
3473 196 092 061 
509618 740 844 
48 964 445 869 
2 991921656 
111984 386 
2 521253 
36 746 
380 



21 155 277 
40 108 145 
32 636 626 
14 939 530 
4 226 027 
762 739 
87 644 
6 238 
261 
6 



J : _ 
285 691 
965 420 
427 969 
821672 
636 156 
720 158 
734 347 
791356 
761786 
025 245 

69 220 
386 



15 954 847 
29 043 492 
22 562 
9 792 558 
2 603 398 
436426 
45833 
2 915 
105 



j =33 
768 088 
273524 
009 720 
415 254 
119525 
275281 
115 884 
994 723 
936 834 
004 377 

16805 
52 



3=34 

11823018 778121 
20 638 880 703 615 
15 283132 698 005 
6 276 719 028 345 
1564 280 125 341 
242 715 583 856 
23 180 300 406 
1 308 416 720 
40 682 440 
617 953 
3 491 



3 =35 

8 605 270 288 428 
14 385 457 021771 
10 136 458 047 648 
3 930 854 880 338 
915 803 054 713 
131021687 166 
11316 367 286 
561964 518 
14 753 318 
174 625 
589 



j=36 

6149017377510 
9 828 892 336 538 
6 577 897 351 181 
2 403 021466 599 
521796 689 085 
68 546 193 842 
5 320 896 693 
230 239072 
5 022 558 
44 548 
73 



j=37 

4 311436 549 811 
6 578 623 394 399 
4173019 367 654 
1432 512 031399 
288 971661588 
34 695 422122 
2 403 511898 
89 612 569 
1 593 264 
10 031 



Appendix p. 11 



j = 38 j=39 j = 40 j = 41 j =42 j =43 j =44 j =45 j = 46 



z yo4 1/4 oyo as / 


1997 470 496 058 


1317 887 697 673 


850 665 829 404 536 654 329 284 


rig fiofi rMi 

J3U 3ZO OZO U41 


198 505 264 692 


116 08' 


5 433 466 66 010 296 070 


4 310 052 558 697 


2 761684 207 395 


1728990 256 047 


1056 510 698 738 629 363 523 349 


365 006 585 603 


205 792 941 637 


112 609 389 066 59 693 833147 


2 585 645 699 396 


1563111427 998 


920894742411 


528 044 208 920 294 272153 267 


159130 489 574 


83 350 267190 


42 202 405 360 20 608 991452 


831790371315 


469 846 262 712 


257 819 091781 


137 219 085 706 


70 712 799 890 


35 214 172 506 


16 908 923 264 


7 809196 055 


3 458 946 792 


155 326 381067 


80 906 285 732 


40 765 822 575 


19 829960 250 


9 291403 735 


4182 750 513 


1803 797 937 


742 649 893 


290 757 073 


16957 718111 


7 985 680153 


3 614180128 


1567 456 613 


649 222 792 


255 780 213 


95 399 665 


33 492 231 


10 990 747 


1039 924 955 


429 469 516 


168 586 789 




62 586 269 


21837 068 


7 105 191 


2134 376 




584 182 


143 089 


32 969 506 


11395 789 


3 672196 




1 092 398 


296 134 


71881 


15 238 




2 717 


383 


466 505 


124 530 




29 815 




6 261 


1117 


161 


17 




1 




1925 


297 




33 




2 














j=47 


j=4S 


j =49 




j = 50 


j=51 


j =52 


j =53 


j =54 


j =55 


j =56 


j =57 


36 432 673 762 


19 480 581012 10 069 527 361 


5019400 819 


2 406133 696 


1 105 659 524 


485 226 772 


202 490 441 79 941549 


29 674 484 


10 280 210 


30 590 692119 


15119155 354 7187 420 893 


3 276 272 681 


1426 882 777 


591 250 737 


231 936 330 


85 622 599 29 531012 


9 430 560 


2 757 016 


9681466 532 


4 362140 777 1878 612 930 


770 201912 


299 178 401 


109 478 461 


37 477 493 


11898 600 


3 465 267 


912 617 


213215 


1464 504 260 


590 430 015 


225 629 513 


81283 454 


27 422 357 


8592 757 


2 474 897 


646 355 


150277 


30 308 


5101 


107 745798 


37 582 223 


12 256 579 


3 706 816 


1029033 


258 783 


57 926 


11258 


1830 


234 


21 


3342 402 


931 895 


234 974 




52 641 


10 232 


1668 


216 


20 


1 






30575 

38 


5 485 
2 


776 




77 


4 














1=58 j = 


:50 j=60 j 


= 61 j = 


62 j 


= 63 j 


= 64 j=65 


j =66 













3293359 964438 254353 59228 11848 1956 250 22 1 

727123 169 643 34 073 5 661 730 65 3 



43020 7 205 940 85 4 

673 62 3 



